u 
o 



A weak spectral condition for the controllability of the 
bilinear Schrodinger equation with application to the 
control of a rotating planar molecule 



U. Boscain abc , M. Caponigro de , T. Chambrion ed , and M. Sigalott 



jC b 



a Centre National de la Recherche Scientifique (CNRS) 
h CMAP, Ecole Poly technique, Route de Saclay, 91128 Palaiseau Cedex, France 
C INRIA, Centre de Recherche Saclay, Team GECO 
d INRIA, Centre de Recherche Nancy - Grand Est, Team CORIDA 
e Institut Elie Cartan, UMR 7502 Nancy-Universite/CNRS, BP 239, Vandceuvre-les- Nancy 54506, 
"q ' France 

o" 

in 

(N 



Abstract 



In this paper we prove an approximate controllability result for the bilinear Schrodinger 
equation. This result requires less restrictive non-resonance hypotheses on the spectrum 
■ of the uncontrolled Schrodinger operator than those present in the literature. The con- 

trol operator is not required to be bounded and we are able to extend the controllability 
result to the density matrices. The proof is based on fine controllability properties of the 
finite dimensional Galerkin approximations and allows to get estimates for the L 1 norm 
of the control. The general controllability result is applied to the problem of controlling 
the rotation of a bipolar rigid molecule confined on a plane by means of two orthogonal 
external fields. 



CO 

1 Introduction 

In this paper we are concerned with the controllability problem for the Schrodinger equation 

^; i d ^ = {H Q + u{t)H l )^. (1.1) 

b : 

Here ip belongs to the Hilbert sphere of a complex Hilbert space H and Hq, Hi are self-adjoint 
operators on %. The control u is scalar- valued and represents the action of an external field. 
The reference model is the one in which H = —A + V(x), Hi = W(x), where x belongs 
to a domain D C R n with suitable boundary conditions and V, W are real-valued functions 
(identified with the corresponding multiplicative operators) characterizing respectively the 
autonomous dynamics and the coupling of the system with the control u. However, equation 
eq-0 can be used to describe more general controlled dynamics. For instance, a quantum 
particle on a Riemannian manifold subject to an external field (in this case A is the Laplace- 
Beltrami operator) or a two-level ion trapped in a harmonic potential (the so-called Eberly 
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and Law model [T5| EE E])- In the last case, as in many others relevant physical situations, 
the operator H Q cannot be written as the sum of a Laplacian plus a potential. 

Equation (11. ip is usually named bilinear Schrddinger equation in the control community, 
the term bilinear referring to the linear dependence with respect to ip and the affine dependence 
with respect to u. (The term linear is reserved for systems of the form x = Ax + Bu(t).) The 
operator Hq is usually called the drift. 

The controllability problem consists in establishing whether, for every pair of states i(jq an d 
tpi, there exist a control it(-) and a time T such that the solution of eq-0 with initial condition 
ip(0) = ipo satisfies ip(T) = Unfortunately the answer to this problem is negative when 
% is infinite dimensional. Indeed, Ball, Marsden, and Slemrod proved in [3] a result which 
implies (see [30]) that equation eq-0 is not controllable in (the Hilbert sphere of) % Moreover, 
they proved that in the case in which H Q is the sum of the Laplacian and a potential in a 
domain D of R n , equation eq-0 is neither controllable in the Hilbert sphere S of L 2 (D, C) nor 
in the natural functional space where the problem is formulated, namely the intersection of S 
with the Sobolev spaces H 2 (D, C) and Hq(D, C). Hence one has to look for weaker control- 
lability properties as, for instance, approximate controllability or controllability between the 
eigenstates of H (which are the most relevant physical states). 

However, in certain cases one can describe quite precisely the set of states that can be 
connected by admissible paths. Indeed in jH [5] the authors prove that, in the case in which 
Hq is the Laplacian on the interval [—1,1], with Dirichlet boundary conditions, and H\ is 
the operator of multiplication by x, the system is exactly controllable near the eigenstates in 
H 7 (D, C) Pi S (with suitable boundary conditions). This result was then refined in [5], where 
the authors proved that the exact controllability holds in H 3 (D, C) D S, for a large class of 
control potentials (see also |23|). 

In dimension larger than one (for Hq equal to the sum of the Laplacian and a potential) 
or for more general situations, the exact description of the reachable set appears to be more 
difficult and at the moment only approximate controllability results are available. 

In [TT] an approximate controllability result for eq-0 was proved via finite dimensional 
geometric control techniques applied to the Galerkin approximations. The main hypothesis 
is that the spectrum of Hq is discrete and without rational resonances, which means that the 
gaps between the eigenvalues of Hq are Q-linearly independent. Another crucial hypothesis 
appearing naturally is that the operator Hi couples all eigenvectors of H . 

The main advantages of that result with respect to those previously known are that: i) it 
does not need H to be of the form —A + V; ii) it can be applied to the case in which Hi is 
an unbounded operator; iii) the control is a bounded function with arbitrarily small bound; 
iv) it allows to prove controllability for density matrices and it can be generalized to prove 
approximate controllability results for a system of Schrddinger equations controlled by the 
same control (see |10|). 

The biggest difficulty in order to apply the results given in [TT] to academic examples 
is that in most of the cases the spectrum is described as a simple numerical series (and 
hence Q-linearly dependent). However, it has been proved that the hypotheses under which 
the approximate controllability results holds are generic [T9| [2U [25]. Notice that writing 
Hq + u(t)Hi = (Hq + eHi) + (u(t) — e)Hi and redefining (u(t) — e) as new control may 
be useful. As a matter of fact, perturbation theory permits often to prove the Q-linearly 
independence of the eigenvalues of (Hq + eHi) for most values of e. This idea was used in |11] 
to prove the approximate controllability of the harmonic oscillator and the 3D potential well 
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Figure 1: Each vertex of the graph represents an eigenstate of i^o (when the spectrum is not 
simple, several nodes may be attached to the same eigenvalue). An edge links two vertices 
if and only if Hi connects the corresponding eigenstates. In this example, (<f>i, Hi<p 2 ) and 
(<j)i, H X (j) 3 ) are not zero, while (0i,#i0 4 ) = (02, #103) = (02,#i04) = 0. 



with a Gaussian control potential. 

Results similar to those presented in [TT] have been obtained, with different techniques, 
in [22] (see also |15[ El I2TJ [23]). They require less restrictive hypotheses on the spectrum of 
Hq (which is still assumed to be discrete) but they do not admit Hi unbounded and do not 
apply to the density matrices. However, it should be noticed that [22J proves approximate 
controllability with respect to some Sobolev norm H s , while the results given in [11] permit to 
get approximate controllability in the weaker norm L 2 . As it happens for the results in |llj . 
the sufficient conditions for controllability obtained in [22] are generic. 

Fewer controllability results are known in the case in which the spectrum of H is not 
discrete. Let us mention the paper |2U| . in which approximate controllability is proved between 
wave functions corresponding to the discrete part of the spectrum (in the ID case), and j 15 j . 

In this paper we prove the approximate controllability of eq-0 under less restrictive hy- 
potheses than those in [11] . More precisely, assume that H has discrete spectrum (Xk)keTSi 
(possibly not simple) and denote by 0^ an eigenvector of H corresponding to in such a way 
that (0fc)fceN is an orthonormal basis of H. Let S be the subset of N 2 given by all (ki, k 2 ) such 
that {4>ku Hi(f>k 2 ) ^ 0. Assume that, for every (j, k) G H such that j ^ k, we have Xj ^ \k 
(that is, degenerate energy levels are not directly coupled by Hi). We prove that the system 
is approximately controllable if there exists a subset S of S such that the graph whose vertices 
are the elements of N and whose edges are the elements of S is connected (see Figure [JJ) and, 
moreover, for every (ji, j 2 ) G S and every (ki,k 2 ) € S different from (ji,j 2 ) and (j 2 ,ji), 

I An — Awl 7^ l-^fci — ^fc 2 l- (1-2) 

As in [TT], Hi is not required to be bounded and we are able to extend the controllability 
result to the density matrices and to simultaneous controllability (see Section 12.21 for precise 
definitions). This extension is interesting in the perspective of getting controllability results 
for open systems. 

Interesting features of our result are that it permits to get L l estimates for the control laws 
and that it does not require the spectrum of A to be simple. Moreover, beside requiring less 
restrictive hypotheses, this new result works better in academic examples where very often 
one has a spectrum which is resonant, but a lot of products (0/^, Hi^k 2 ) which vanish. The 
consequence of the presence of these vanishing elements (i.e., of the smallness of S) is that 
less conditions of the type (jl.2p need being verified. 

The condition on the spectrum given above is still generic and it is less restrictive than 
the one given in [22], which corresponds to the case S = {(k , k) \ k G N, k 7^ k } for some 
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Figure 2: The bipolar rigid molecule confined to a plane. 



k G N and where condition (ll.2p is required for every (k\, k 2 ) G N 2 \{(j!, j 2 ), ji)}- Notice 
however that the approximate controllability result given in [22J is still in a stronger norm. 

The idea of the proof is the following. We recover approximate controllability for the 
system defined on an infinite dimensional Hilbert space through fine controllability properties 
of the N- dimensional Galerkin approximations, N G N, which allow us to pass to the limit 
as N — > oo. More precisely, we prove that, for n, N G N with N ^> n ^> 1, for given initial 
and final conditions ipQ, ipi in the Hilbert sphere of % which are linear combinations of the 
first n eigenvectors of Hq, it is possible to steer ipo to ipi in the Galerkin approximation of 
order N in such a way that the projection on the components n + 1, . . . , N has arbitrarily 
small norm along the trajectory. This kind of controllability for the Galerkin approximation 
of order N is proved in two steps: firstly, thanks to a time-dependent change of variables we 
transform the system in a driftless one, nonlinear in the control, and we prove the result up 
to phases. The change of variables was already introduced in [TJ [TT|; the technical novelty of 
this paper is the convexification analysis for the transformed system, which allows to conclude 
the controllability with less restrictive non-resonance hypotheses. Secondly, the control of 
phases is obtained via a classical method, using as pivot an eigenstate of H and exploiting 
the controllability (up to phases) of the time-reversed Schrodinger equation. This last step 
requires some further arguments in the case of simultaneous controllability. 

In the second part of the paper we apply our result to the problem of controlling a bipolar 
rigid molecule confined on a plane by means of two electric fields constant in space and 
controlled in time, oriented along two orthogonal directions (see Figure [2]). The corresponding 
Schrodinger equation can be written as 

l ^T L= (-^ + Mt)™*(0) + Mt)M0))w,t), 9e$\ (1.3) 
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where S 1 = R/27rZ. Notice that a controlled rotating molecule is the most relevant physical 
application for which the spectrum of H is discrete. 

The system described by (11 .3p is not controllable if we fix one control to zero. Indeed by 
parity reasons the potential cos (9) does not couple an odd wave function with an even one 
and the potential sin (9) does not couple wave functions with the same parity. 

Up to our knowledge, the controllability result presented in this paper is the only one 
which can be directly applied to system eq-1. Indeed, both the results obtained in jTT] and |22] 
seem to require sophisticated perturbation arguments in order to conclude the approximate 
controllability of eq-1. 

The proof of the controllability of eq-1 by means of the controllability result obtained in 
this paper is not immediate since we have to use in a suitable way the two controls. The idea 
is to prove that, given an initial condition ip Q which is even with respect to some 9 G S 1 , by 
varying the control u = (wi,^) € R 2 along the line R(cos(#), sin(#)), it is possible to steer it 
(approximately) towards any other wave function even with respect to 9. In particular, it is 
possible to steer any eigenfunction (which is necessarily even with respect to some 9) to the 
ground state (that is, the constant function 1/a/27t), which can, in turn, be steered towards 
any other eigenfunction. The argument can be refined to prove approximate controllability 
among any pair of wave functions on the Hilbert sphere. 

The structure of the paper is the following. In Section [2] we introduce the class of systems 
under consideration and we discuss their well-posedness. Then, we state the main results 
contained in the paper. Section [3] is devoted to the case in which % is finite dimensional. 
Sections SJ \5[ and contain the proof of the main results and in Section 15.61 we present 
estimates on the L 1 norm of the control. Section [7] contains an application to the infinite 
potential well, showing controllability and establishing L 1 estimates of the control. Section [8] 
provides the application to the bipolar planar molecule evolving on the plane. 



2 Framework and main results 
2.1 Settings and notations 

As in [TT], we use an abstract framework instead of a presentation in terms of partial differ- 
ential equations. The advantage of this presentation is that it is very versatile and applies 
without modification for Schrodinger equation on a (possibly unbounded) domain of R n or on 
a manifold such as S 1 (see Section [8]). To avoid confusion, let us stress that H Q and Hi, intro- 
duced in the introduction, are self-adjoint operators while A = —iH and B = —iH\, used in 
what follows, are skew-adjoint. Hereafter N denotes the set of strictly positive integers. We 
also denote by U("H) the space of unitary operators on TL. 

Definition 2.1. Let % be an Hilbert space with scalar product (•, •} and A, B be two (possibly 
unbounded) linear operators on with domains D(A) and D(B). Let U be a subset of R. 
Let us introduce the formal controlled equation 

^(t) = (A + u(t)B)4>(t), u(t) e U. (2.1) 
We say that (A, B, U, $) satisfies (21) if the following assumptions are verified: 
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(211) $ = ((pkjke'N is an Hilbert basis of % made of eigenvectors of A associated with 
the family of eigenvalues (i\k)k£~N] 

(212) (p k G D(B) for every fceN; 

(213) A + uB : span{0fc | k G N} — > H, is essentially skew-adjoint for every u £ U; 

(214) if j ^ k and Aj = X k then £0 fc ) = 0. 

Remark 2.2. If A has simple spectrum then (214) is verified. If all the eigenvalues of A 
have finite multiplicity, then, up to a change of basis, hypothesis (214) is a consequence of 
(211-2-3). 

A crucial consequence of assumption (213) is that, for every constant u in. U, A + uB 
generates a group of unitary transformations e t (. A+uB ) : % — >. The unit sphere of "H is 
invariant for all these transformations. 

Definition 2.3. Let (A, B, U, $) satisfy (21) and u : [0,T] — > U be piecewise constant. The 
solution of 02. ip with initial condition tp G H is 

iW) = TWo), (2.2) 

where T u : [0,T] — > U(H) is the propagator of ( 12. ip that associates, with every t in [0,T], the 
unitary linear transformation 

T u = e^-Si'')^ 5 ) o e *i-i(^+ u i-i B ) o • • • o ^(A+u.B)^ 

where U<t < Y?i=i *i and n ( r ) = u i if ^ - r < SLi 

The notion of solution introduced above makes sense in very degenerate situations and can 
be enhanced when B is bounded (see [3] and references therein). 
Note that, since 

for every n G N, ipo £ H, and u £ U, then, for every solution ip(-) of ( 12. ip . the function 
t H- (i[>(t),(j) n ) is absolutely continuous and satisfies, for almost every t G [0, T], 

j t V(O) = - ({A + u{t)B)<f> n , m) ■ (2-3) 

2.2 Main results 

As already recalled in the introduction, exact controllability is hopeless in general. Several 
relevant definitions of approximate controllability are available. The first one is the standard 
approximate controllability. 

Definition 2.4. Let (A, B,U, $) satisfy (21). We say that (12. ip is approximately controllable 
if for every ipo,ij)i in the unit sphere of "H and every e > there exist a piecewise constant 
control function u : [0,T] — > U such that — T^(ijjo)\\ < e. 

Recall that A has purely imaginary eigenvalues (i\k)k&N with associated eigenfunctions 
(</>fc)fceN- Next we introduce the notion of connectedness chain, whose existence is crucial for 
our result. 
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Definition 2.5. Let (A, B, U, $) satisfy (21). A subset S of N 2 couples two levels j, k in N, 
if there exists a finite sequence ((s\, s^), . . . , (s?, S2)) in S* such that 



s\ = j and S2 



~3 • 



for every 1 < j < p — 1 ; 
i) ^ for 1 < j <p. 



S is called a connectedness chain (respectively m-connectedness chain) for (A, B, U, $) if 
S (respectively 5 PI {1, ... , m} 2 ) couples every pair of levels in N (respectively in {1, . . . , m}). 

A connectedness chain is said to be non-resonant if for every (si,s 2 ) in S, \X Sl — X S2 \ 7^ 
|A tl - At 2 | for every (ti,t 2 ) in N 2 \ {(si,s 2 ), (s 2 ,si)} such that (4> t2 , B(j) tl ) ^ 0. 

Theorem 2.6. Let 5 > and let (A,B, [0,5],$) satisfy (21). // i/iere exists a non-resonant 
connectedness chain for (A,B, [0,5],$) then (12. ip is approximately controllable. 



Theorem 12.61 is a particular case of Theorem 12. 11^ stated in the next section. 

Remark 2.7. Notice that in the assumptions of Theorem 12.61 we do not require that the 
eigenvalues of A are simple. Take for instance T-L = C 4 , U = [0, 1], and 
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A connectedness chain is given by {(1, 2), (2, 1), (1, 3), (3, 1), (2, 4), (4, 2)}. The corresponding 
eigenvalue gaps are |A 2 — Ai| = 1, | A3 — Ai| = 3, and | A4 — A 2 | = 2. Hence, the connectedness 
chain is non- resonant . 

The following proposition gives an estimate of the L 1 norm of the control steering (12. ip 
from one eigenvector to an e-neighborhood of another. A generalization of this proposition is 
given by Theorem 12.131 

Proposition 2.8. Let 5 > 0. Let (A,B, [0,5],$) satisfy (21) and admit a non-resonant chain 
of connectedness S. Then for every e > and (j,k) G S there exist a piecewise constant 
control u : [0,TJ — > [0,5] and 9 G R such that ||Y^(0j) — e* 9 0fc|| < e and 

5ir 



2.3 



Simultaneous controllability and controllability in the sense of 
density matrices 



We define now a notion of controllability in the sense of density matrices. Recall that a density 
matrix p is a non-negative, self-adjoint operator of trace class whose trace is normalized to 
one. Its time evolution is determined by 

p(t) = T t >(0)Tf 

where T"* is the adjoint of T". Notice that the spectrum of p(t) is constant along the motion, 
since, for every t, p(t) is unitarily equivalent to p(0). 
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Definition 2.9. Let (A, B, U, $) satisfy (21). We say that (12. ip is approximately controllable 
in the sense of the density matrices if for every pair of unitarily equivalent density matrices 
po, pi and every e > there exists a piecewise constant control u : [0, T] — > U such that 

|| Pl -T>oTT|| <e, 

in the sense of the operator norm induced by the Hilbert norm of H. 

Definition 2.10. Let (A, B, U, $) satisfy (21). We say that (12. ip is approximately simultane- 
ously controllable if for every r in N, . . . , ip r in %, T in U("H), and e > there exists a 
piecewise constant control u : [0, T] — > U such that, for every 1 < k < r, 



< e. 



The following result is proved in Sections HI [5J and [6j 

Theorem 2.11. Let 5 > and let (A,B, [0,5],$) satisfy (21). If there exists a non-resonant 
connectedness chain for (A, B, [0, S], $) ; i/ien (12. ip is approximately simultaneously control- 
lable. 

Simultaneous controllability implies controllability in the sense of density matrices (see 
Proposition lA.ip . Hence we have the following. 

Corollary 2.12. Let 5 > and let (A,B, [0,5],$) satisfy (21). If there exists a non-resonant 
connectedness chain for (A,B, [0,5],$), then (12. ip is approximately controllable in the sense 
of the density matrices. 

Theorem 2.13. Let (A, B, U, $) satisfy (21) and admit a non-resonant chain of connectedness. 
Then there exists a basis $ = (0/«)feeN of eigenvectors of A and a subset S of N 2 such that, 
for every m G N, S is a m- connectedness chain for (A, B,U,&) . Moreover, let 5 > and 
U = [0,5], then for every e > and for every permutation a : {l,...,m} — > {!,..., m} 
there exist a piecewise constant control u : [0, T u ] — > [0, 5] and 8\, . . . , 9 m in R for which the 
propagator T" of (12. ip satisfies ||Y^ u <^ — e l0 '<j) a (i) \\ < e for every 1 < I < m and 

5 n (2 m ^ 1 - 1) 

MU 1 < 



AM{\U k ,B^)\ : (j,k)eS, l<j,k<m} 



The proof of the first part of the statement of Theorem 12.131 is given in Section 14.41 The 
second part is proved in Section 15.61 

Notice that a lower bound on the L 1 norm of the control was already proved in (see 
Proposition I5.10p . 



3 Finite dimensional case 

Denote by u(n) and su(n) the Lie algebras of the group of unitary matrices U(n) and its 
special subgroup SU(n) = {M e U{n) \ det M = 1} respectively. 

Here we address the case where "H is of finite dimension n. Equation (12. ip then defines 
a bilinear control system on U(n). Finite dimensional systems of the type (12.11) have been 
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extensively studied. A necessary and sufficient condition for controllability on SU(n) (i.e. the 
property that every two points of SU (n) can be joined by a trajectory in U(n) of system (12. ip ) 
is that the Lie algebra generated by A and B contains su(n). This criterion is optimal, yet 
sometimes too complicated to be checked for n large. Easily verifiable sufficient conditions for 
controllability on SU(n) have been thoroughly studied in the literature (see for instance |13| 
and references therein). Next proposition gives a new sufficient condition, slightly improving 
those in |3U) and [TTj Proposition 4.1]. Its proof is based on the techniques that we extend to 
the infinite dimensional case in the following sections. 

The controllability result is obtained under a slightly weaker assumption than (21). 

Proposition 3.1. Let H = C n . Let (A, B,U,&) satisfy (211 — 2 — 3) and admit a non- 
resonant connectedness chain S. Assume, moreover, that Xj ^ X k for every (j,k) G S. Then 
the control system (12.1 p is controllable both on the unit sphere of C n and on SU(n), provided 
that U contains at least two points. If, moreover, tr A ^ or tr B ^ 0, then the control system 
(12. ip is controllable on U(n). 

Proof. For every 1 < j, k < n, let be the n x n matrix whose entries are all zero, but the 
one at line j and column k which is equal to 1. We denote by a^ k and bj k the (j, k)-th entry 
of A and B, respectively. 

Recall that, for any two n x n matrices X and Y, &dx{Y) = [X, Y] = XY — YX, and 
compute the iterated matrix commutator 

n 

ad^(#) = ( a ii - a kk) p bjke$ ■ 
j,k=i 

Fix (j,k) in S. By hypothesis, for every l,m in {l,...,n} such that {l,m} ^ {j, k}, 
) 2 7^ ( a u - a-mm) 2 or bi m = 0. There exists some polynomial P jk with real coefficients 
such that Pjk{{ajj - a kk ) 2 ) = 1 and P jk ({a u - a mm ) 2 ) = if [a^ - a kk ) 2 ^ (an - a mm ) 2 . Let 
Pjk = Yi=oChX h . Then 

J2c h ^ 2 A h (B)=b jk e^ + b kj e^. 

h=0 

As a consequence, 

d 

Ch^ 2 A +1 {B) = (ajj - a kk ) (b jk ef k ] + 6^4?) = *( A j ~ X k) (fyk^ + b jk e$) , 



h=0 



and then the two elementary Hermitian matrices e$ — ejy and ie& + iejy also belong to 
Lie(y4, B). Because of the connectedness of B and thanks to the relation 



>) >) 



one deduces that su(n) C Lie(A,B). 

If tr A = ti B = 0, then A and B belong to su(n), hence su(n) = Lie(A, B). If tr A ^ 
or tr B 7^ 0, then A or B does not belong to su(n) and u(n) = Lie(A, B). This completes the 
proof of the controllability of the control system (12. ip on SU(n) and U(n). 
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It remains to prove the controllability on the unit sphere S n of C n . Fix Xq,Xi in S n , and 
consider an element of g x G SU(n) such that giXo = X\. According to what precedes there 
exists a trajectory g in U(n) of ( 12. ip from I n to g\. The curve t t— > g(t)xo is a trajectory 
of (EH) in S n that links x to x\. □ 



4 Convexification procedure 

Sections HI |5l and |6] are devoted to the proof of Theorem 12.111 in the case in which Ti has 
infinite dimension. 



4.1 Time-reparametrization 

We denote by PC the set of piecewise constant functions u : [0, oo) — > [0, oo) such that there 
exist u±, . . . , Up > and — t± < ■ ■ ■ < t p+ \ = T u for which 

p 

u:t^r ^ w iXfe,t i+ i)(*)- 

Let us identify u = Y^j=i u jX[t Jt t J+1 ) with the finite sequence (uj, Tj)i<j<p where Tj = tj + \ — tj 
for every 1 < j < p. 
We define the map 

V: PC ^ PC 

which satisfies the following easily verifiable properties. 

Proposition 4.1. For every u E PC, V o V(u) = u and \\V(u)\\ L i = Yl P i=i T j- 

Assume that (A, B, U, $) satisfies (21). In analogy with Definition 12. 3[ we define, for every 
u = Y^=i u jX[tj,tj +1 ) £ PC such that u{t) G U for every t > 0, the solution of 

^(t) = (u(t)A + B)m, (4-1) 

with initial condition tpo G H as 

^(t) = gCt-tOtaA+B) Q . . . Q e *i(«i^+s)(^ ) s 

where U <t < ti + i. 

System (14.11) is the time reparametrization of system (12. ip induced by the transformation 
V, as stated in the following proposition. 

Proposition 4.2. Let u = {uj,Tj)i<j< p belong to PC and ipo be a point ofH. Let ip be the 
solution of Ii2.1\) with control u and initial condition tp 0; and if; be the solution of jjj4-l\ ) with 
control V(u) and initial condition ipQ. Then ip (T u ) = ip (||m||x,i). 

Proof. It is enough to remark that, if u ^ 0, for every t G [0, oo), e *( A+uB ) = e tu {u A+B ) . □ 

As a consequence of Proposition 14.21 it is equivalent to prove controllability for (12. ip with 
U = (0, 5] or to prove controllability for system (14. ip with control u G [1/5, oo). 
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4.2 Convexification 

For every positive integer N let the matrices 

A™ = diag^Ax, . . . , t\ N ) and B« = £ty fc »£ fc=1 =: (b jk )f )k=1 , 

be the Galerkin approximations at order N of A and 5, respectively. Let t h- >• ^(t) be a 
solution of 

corresponding to a control function u and consider t>(t) = J^u(r)dT. Denote by d(B) the 
diagonal of flW and let flW = flW - d(5). Then g : t ^ e~ v ^ A(N) - td ^ip(t), is a solution 
of 

q(t) = e -v(t)AW-t d (B)£ { N) e v(t)AW + td(B) q ( t) _ (9jv) 

Let us set 

= e -«AW-td(B)^(JV) e «AW+td(fl)_ (42) 

Lemma 4.3. Let K be a positive integer and 71, ... , 7^ G R\ {0} 6e sitc/i £/iat |7i| 7^ |7j| /or 
j = 2, . . . , K. Let 

¥>(f) = (e*» ....e**). 

Then, for every to G R, we /iai>e 



conv(^([t ,oo)) D z/S 1 x {(0, ...,0)}, 

where v = YIT=2 cos (If) > ^- Moreover, for every R > and £ G S 1 i/iere exzsfo a sequence 
(^fc)feGN sitc/i £/iat —tf c >R and 



1 - 

lim - V^) = K,0,...,0). 
h-s-oo h z — ' 
k=l 



Proof. Since 

tp{t - t ) = ( e - it07l e ii71 , . . . , e -^o7K e ^) ) (4.3) 

it is enough to prove the lemma for t = 0. We can suppose that \ jx\ = 1 and, up to a reordering 
of the indexes, that there exist n and n such that 1 < n < n < K, ^ for every 
i,j G {1, . . . ,n}, 72, • • • ,7n e Z, 7n + i, ... ,7^ G R\Z, and {|7„+i|, • • • , |7n|} C {I72I, • • •, |7n|}- 
Consider the 2 n_1 real numbers defined as follows: let 

h = 0, 

and for k G {1, . . . , n — 1} and j G {1, . . . , 2 fc ~ 1 }, 



l7fc+i| 

Up to a reordering of the tj, we can suppose that = t\ < fa < ■ ■ ■ < fan-i. Take an integer 
r larger than R/2tt, then set tj = tj + 2nr(j — 1), in such a way that t^ — tfc_i > R for every 
A; = 2,...,2 n ~ 1 . 
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Now consider the arithmetic mean of the l-th (complex) coordinates of <p(ti), . . . , (p(t 2 n-i). 
We show that this quantity is zero for I = 2, . . . , h. Indeed, from the definition of tj, we have 

2™- 1 n-1 

E e^ 7i = Y[ (1 + e ^ 7i/l7fc+l1 ) , 

j=l k=l 

which is zero since so is the fc-th factor when |7/| = |7fc+i|. 

On the other hand, the arithmetic mean of the first coordinate is uniformly bounded away 
from zero. Indeed 



- exp f £ log ( cos (It) ) ) = 27 • ( 4 - 4 ) 

Since log (cos (£)) ~ — ^3 as k tends to infinity, then the sum ^2 k>2 ^°S (cos (^)) converges 
to a (negative) finite value /. As a consequence, v = exp(Z) is a positive number. 

Therefore we have found a sequence of numbers tj such that the arithmetic mean of the 
first coordinate of </?(ii), . . . , (p(t 2 n-i) is uniformly bounded away from zero and the arithmetic 
means of following h — 1 coordinates are zero. According to ( 14. 3p . the role of ti, . . . , £2™-* can 
equivalently be played, for every k G N, by the 2 n ~ 1 -uple 

tj = tj + 27rmA; , j = 1, . . . , 2 n_1 , 

where the integer m is larger than r + t 2 ™-i /27r. Now, let I 6 {n + 1, . . . , i^}, so that 7; ^ Z. 
For every /isN, the arithmetic mean of the Z-th coordinate of the points (p(tj) (k — 0, . . . , h, 
j = l,...,2 n ~ 1 ) is 

2"- 1 h 2"- 1 h 

2"" 1 (/z + l)^^ 2 n - 1 (h+l)^ 

v ' j=l k=0 v ' j=l fc=0 

1 ^» \ 1 1 — p i27rm(/i+l)T; 



3=1 



+ 1-e 



Therefore, we found a sequence of points in the convex hull of </?([0,oo)) converging to 
(2 1_n 5Ty=i e* 71 ^ , 0, . . . , 0). The lemma follows from (14. 4p and by rotation invariance (see ( 14.3 j) ) . 



□ 



2 4 

Remark 4.4. In order to estimate u, notice that, for every x in (—1, 1), — ^ — < log(cos(x)). 
Hence, taking x — for k >2, 



00 ^ 4 00 



7T 4 1 vr 8 + 240 vr 4 - 1980 vr 2 



E (cos (^) ) > -y E ^ - E F 15840 

from which one deduces v > exp (— 7r8+240 1 gg~ 1980 71-2 j > |. Numerically, one finds v ~ 0.430. 
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4.3 An auxiliary system 

Let (A, B, U, $) satisfy (21). With every non-resonant connectedness chain S for (A, B, U, <3>) 
and every n G N we associate the subset 

S n = {(j,k)eS\l<j,k<n,j^k} 

of S and the control system on C n 

where 9 = 9(t) 6 S 1 and (j, k) = k(t)) G are piecewise constant controls. Recall that 

is the n x n matrix whose entries are all zero but the one of index (j, k) which is equal to 

1 and that v = Yl'k'=2 cos (ik) ( see Lemma f4 . 3 1) . 

The control system (E n ) is linear in x. For every 9 in S 1 and every 1 < j, k < n, j ^ k, the 
matrix e t9 e^ — &~ %6 &\^ is skew-adjoint with zero trace. Hence the control system (E n ) leaves 
the unit sphere S n of C n invariant. In order to take advantage of the rich Lie group structure 
of group of matrices, it is also possible to lift this system in the group SU(n), considering x 
as a matrix. 

4.4 Existence of a n-connectedness chain 

Notice that system (S n ) cannot be controllable if 5* is not a n-connectedness chain (see |114 
Remark 4.2]). This motivates the following proposition. 

Proposition 4.5. Let (A, B,U,Q) satisfy (21). // there exists a connectedness chain for 
(A,B,U,Q), then there exists a bijection a : N — > N such that, setting $ = (4> a (k))keN 
and S = {(a(j),a(k)) : (j,k) G S}, (A, B,U,&) satisfies (21) and S is a n-connectedness 
chain for (A, B,U,&) for every n mN. 

Proof. Following |19[ Proof of Theorem 4.2], a can be constructed recursively by setting 
<j(l) = 1 and a(n + 1) = mina ({<r(l), . . . ,a(n)}), where, for every subset J of N, a(J) = 
{k G N \ J | B(f) k ) ^ for some j in J} . □ 

Proposition 14.51 proves the first part of the statement of Theorem 12.131 

5 Modulus tracking 

The aim of this section is to prove the following proposition, which is the main step in the 
proof of Theorem 12.61 

Proposition 5.1. Let 5 > 0. Let (A,B, [0,5],$) satisfy (21) and admit a non-resonant con- 
nectedness chain. Then, for every continuous curve T : [0, T] — > U("H) such that Tq = 1%, 
r in N, and e > 0, there exist T u > 0, a continuous increasing bijection s : [0,T] — > [0,T U ], 
and a piecewise constant control function u : [0,T U ] — > [0,5] such that the propagator T u of 
equation (12. ip satisfies 

iK^t^l-K^T^H <e for every t G [0,T], j G N, 

for every G span{0 l5 . . . , r } with ||0|| = 1. 
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The proof of Proposition 15.11 splits in several steps. In Section 15.11 we recall some classical 
results of finite dimensional control theory, which, in Section 15.2} are applied to system (E n ) 
introduced in Section fl~3l In Section I5TB1 we prove that system (Ojv) can track in projection 
the trajectories of system (S n ). Then we prove tracking for the original infinite dimensional 
system in Section 15.41 The proof of Proposition 15.11 is completed in Section 15.51 

5.1 Tracking: definitions and general facts 

Let M be a smooth manifold, U be a subset of R, and / : M x U — >■ TM be such that, for 
every x in M and every u in U, f(x,u) belongs to T X M and f{-,u) is smooth. Consider the 
control system 

x = f(x,u), (5.1) 

whose admissible controls are piecewise constant functions u : R — >■ U. For a fixed u in U, we 
denote by f u the vector field x \— > f(x, u). 

Definition 5.2 (Tracking). Given a continuous curve c : [0, T] — > M we say that system (15.11) 
can track up to time reparametrization the curve c if for every e > there exist T u > 0, an 
increasing bijection s : [0,T] — > [0,T U ], and a piecewise constant control u : [0,T U ] — > U such 
that the solution x : [0, T u ] — > M of ( 15. ip with control u and initial condition x(0) = c(0) 
satisfies dist(x(s(t)), c(t)) < e for every t G [0, T], where dist(-, •) is a fixed distance compatible 
with the topology of M. If s can be chosen to be the identity then we say that system (15.11) 
can track c without time reparametrization. 

Notice that this definition is independent of the choice of the distance dist(-,-). Next 
proposition gives well-known sufficient conditions for tracking. It is a simple consequence of 
small-time local controllability (see for instance [T2"| Proposition 4.3] and |16j). 

Proposition 5.3. If, for every x in M, {f(x,u) \ u G U} = {— f(x,u) \ u G U} and 
Lie x ({f u | u G U}) = T X M, then system 115. 1\) can track up to time reparametrization any 
continuous curve in M. 

5.2 Tracking in (E n ) 

We now proceed with the first step of the proof of Proposition 15.11 Using Proposition 15.31 we 
can prove the following. 

Proposition 5.4. Let (A, B,U,<&) satisfy (21). Let S be a non-resonant connectedness chain 
for (A, B, U, $) such that, for every n in N , S is a n- connectedness chain. Then, for every n 
in N ; the finite dimensional control system (£„) can track up to time-reparametrization any 
curve in SUin). 

Proof. Recall that S n = {(j, k) G S \ 1 < j, k < n,j ^ k}. In order to apply Proposition 15.31 
we notice that the set 

V(s) = {u\b jk \ (e ie e^ - e" ie eg } ) x : 6 G S 1 , (j, k) G S n ] 

is symmetric with respect to and we are left to prove that the Lie algebra generated by the 
linear vector fields x i— > v\bjk\(e ej^ — e~ %e e^j)x contains the whole tangent space 5u(n)x of 
the state manifold SU(n). The latter condition is verified if and only if B^ is connected, as 
shown in the proof of Proposition 13.11 □ 
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5.3 Tracking trajectories of (S n ) in (On) 

Next proposition states that, for every N > n, system (On), defined in Section HI can track 
without time reparametrization, in projection on the first n components, every trajectory of 
system (£„). 

Hereafter we denote by the projection mapping a N x N complex matrix to the n x n 
matrix obtained by removing the last N — n columns and the last N — n rows. 

Proposition 5.5. Let 5 > 0. Let (A, B, [0, 5], satisfy (21) and admit a non-resonant 
connectedness chain S . For every n, N G N, N > n, e > and for every trajectory 
x : [0, T] — > SU (n) of system (E n ) with initial condition x(0) = I n there exists a piece- 
wise constant control u : [0,T] — > [1/5, +oo) such that the solution y : [0, T] — > SU(N) of 
system (Oat) with initial condition In satisfies 



\\x(t) - Ilf \y{t))\\ < e for every t G [0,T]. 

Proof. Given a trajectory x(t) of system (S n ) with initial condition x(0) = /„, denote by 
(j,k) = (j(t),k(t)) G S n and 9 = 9(t) G S 1 its corresponding control functions. Being these 
functions piecewise constant, it is possible to write [0,T] = Up=o[^P'Vfi] ^ n sucn a wa Y that 
j, k, and 9 are constant on [t p , t p+ i) for every p — 0, . . . , q. 

We are going to construct the control u by applying recursively Lemma 14.31 Let 5 > 1/5. 
Fix p G {0, . . . , q} and j, k, 9 such that (j(£), k{t)) = (j, k) and 9{t) = 9 on [t p , t p+ i). Apply 
Lemma S3] with 71 = Xj - X k , {72, . . . ,j K } = {Xi - X m \ l,m G {1, . . . ,N},bi m ^ 0, {l,m} ^ 

-T + 5T, and to = t (p) to be fixed later depending 



{j, k}, and I ^ m}, R = max (ijm)e5n 

on p. Then, for every 77 > 0, there exist h = h(p) > 1/rj and a sequence (t/;^ =1 such that 



I A; — A m | 



w{ > t , w 



P _ 
a 



W' 



.j > R, and such that 



1 h 



Aj)«)S 



if) 

\ b jk\ 



and 



5E 



,i(A r A m )ioS 



< V, 



for every /, m G {1, . . . , N} such that bi m 7^ 0, {I, m} 7^ {j, fc}, and / 7^ m. 

Set t~p = t p + (t p+ i — t p )a/h, a = 0, . . . , h, and define the piecewise constant function 



*>(*) = EE («* 

p=0 a=l ^ 



pbjj-hk 



a 1 a \ \ 
Aj — Ak 



(5.2) 



Note that by choosing to(p) = + -R for p = 1, . . . , q and io(0) = R we have that f^(t) 

is non-decreasing. 

Set M(t) = v\bj( t )k(t) I ( e te ^e[^2i,.u\ — e^^e^L^ I • From the construction of iv, we have 



j(t)k{t) 
t 



-mm ) ■ 



'd N (s,v r) (s))ds 



T?— >0 







M(s)ds, 



(5.3) 
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arctan(5) 



Figure 3: The piecewise constant function v v (in bold) and the piecewise linear approximation 
v„ with slope greater than 8. 



uniformly with respect to t G [0,T], where i9n is defined as in (14. 2ft . This convergence 
guarantees (see for example |2] Lemma 8.2]) that, denoting by y v (t) the solution of system (On) 
with control v v and initial condition In, U ( h f \y v (t)) converges to x(t) as rj tends to uniformly 
with respect to t G [0, T\. Hence, for every rj sufficiently small, 

||nf) o y v (t) - x(t)\\ < £ - for every t G [0,T]. 



If the functions v v were of the type t h-> f*u(s)ds for some u : [0, T] — > [1/5, oo) piecewise 
constant, then we would be done. For every 77 > consider the piecewise linear continuous 
function v v uniquely defined on every interval [t p ,t p+ i) by 



v v (t p ) 
v v (t) 



W 



p-i 



v v (t) 



iHGULi[^-i>d + 

for a = 1, . . . , h, 

\h r_P 1 (tp+v—tp) 



s iHGu: =1 [Ci+ 



where we set = (see Figure [573]) . On each interval [t[_ ± + {lp+ li lp> , if) the difference 
between v v and is bounded in absolute value by 5(t p+ i — t p )/h. Therefore, 



SUp < \\$ N {t,V v (t)) -^N^V^t)) || I t G 



u 




tends to zero as 77 tends to 0. 

Since ||$;y(i, v)\\ is uniformly bounded with respect to (t,v) G [0,T] x R and the measure 
°f U/ P [ r z-i5 r z-i + ^ P+ fe2~ ty ) g° es to as r] goes to 0, we have 
-t 



{■&n(t, v v (r)) - $ n (t, Vt,(t))) dr 



77^0 



uniformly with respect to t G [0,T] . 
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In particular, for rj sufficiently small, if y v denotes the solution of system (Gat) with control 
v v and initial condition 1^, then 

\\y v (t) - y n (t)\\ < | for every t G [0,T]. 

Finally, u can be taken as the derivative of v v , which is defined almost everywhere. □ 

5.4 Tracking trajectories of (E n ) in the original system 

Next proposition extends the tracking property obtained in the previous section from the 
system (©at) to the infinite dimensional system (14.11) . We denote by 11^ : Ti — > C n the 
projection mapping ip G % to ((0i, ip), . . . , ((f) n , ip}) G C n and we write </>j^ for U n (j)k- 

Proposition 5.6. Let 5 > and let (A, B, [0, 5], $) satisfy (21). For every e > 0, n G N, and 

for every trajectory x : [0, T] — > SU (n) of system (£„) with initial condition I n there exists 
a piecewise constant function u : [0, T] —> [1/5, +oo) such that the propagator T" of 04. ip 

S(2^2s/ie5 

||(4 n) ,x(on n 0)|-|(^,Tr(0)>||< £ 

/or every <fi G span{0i, . . . , n } with ||0|| = 1 and every t m [0, T] and j in N. 

Proof. Consider \x > 0. For every j G N the hypothesis that <f>j belongs to D(B) implies 
that the sequence (bjk)k&x is in I 2 . It is therefore possible to choose N > n such that 
Sfc>7v \bjk\ 2 < ^ for every j = 1, . . . ,n. By Proposition I5.5[ for every 77 > and for every 
trajectory x : [0, T] — >■ SU(n) of system (E n ) with initial condition J re , there exists a piecewise 
constant control m 77 : [0, T] — > [1/5, +00) such that the solution y v of system (©at) with initial 
condition ijy satisfies 

||a;(t)-nS^(y'(t))|| <»/. 

Denote by R v (t,s) : C N — > 0^,0 < s,t < T, the resolvent of system (©at) associated 
with the control r T? (t) = J Q * u v (T)dr, so that ^(t) = ^(t, 0). Fix G span{0i, . . . with 
= 1 and set 

QV(t) = e -vO(t)A-td(B) T ui 

The components of Q v {t), say 

g ;(t) =e -^^W-^(^,Tr(0)), jGN, 
satisfy, for almost every t G [0,T], 

00 

gj(t) = J2 bjke^-^V+^-^qlit). (5.4) 
fc=i 

Therefore Q v N (t) = U N Q v (t) = (qi(t), . . . , q r) N {t)) T satisfies the time-dependent linear equation 

Qi(t) = Mt,v v (t))Qi(t)+pm, 

where P v N {t) = (£ fc>Ar b lk e< x ^> n V^™'^ q v k , . . . , J2 k>N W i(Afe " AivKW+ * (6 "- f ' JViv) gD T - 
Hence t 

Q v N (t) = R r '(t,0)U N <f ) + [ R^(s,t)P^s)ds. 

Jo 
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Consider the projection of the equality above on the first n coordinates. Notice that, because 
of the choice of N, the norm of the first n components of P^(t) is smaller than y/Jm. By (15. 3p . 
R v (s, t) converges uniformly, as 77 tends to 0, to a time-dependent operator from into itself 
which preserves the norm of the first n components. Then, there exists 77 sufficiently small 
such that 



R v (s,t)P%(s)ds 



< 2TJ]m. 



Hence 



|iW(t) -*(t)n„0|| < \\n n Q*{t) -nf)(^(t,o))n n 0|| + \\Y[^\w(t^))ii n <t>-x(t)Ti n <p\\ 

< 2T v fjIn~ + ri < J, (5.5) 



if fi < e 2 / (32nT 2 ) and 77 < T^fjm. In particular, for j — 1, . . . , n, 

||(4 n) ,x(t)n^>|-|(^,Tf(0))||<|. 

It remains to prove the statement for j > n. From (15. 5p it follows 

n 9 

Ei"?wi 2 >(!-!) . 

3=1 

then, since T" is a unitary operator for every t, we have 



£l<#)l 2 <i-(i-f) 2 <- 



□ 



5.5 Proof of modulus tracking 

The following proposition allows to reduce the tracking problem stated in Proposition 15.11 to 
the tracking of a curve in SU (n) . 

Proposition 5.7. For every continuous curve T : [0,T] — > U^H), e > 0, and x G N, 

there exist n > x and a continuous curve F n : [0, T] — > SU(n) such that \{<pj 1 T t (f)f t ) — 
(<f>f\ F n (t)0^ } )| < s for every t in [0, T], 1 < k < x, and j G N. 

Proof. For every n in N, define the function g n : t 1— > Ylk=i 2~2™=i \ { ( t ) h ^t<Pk)\ 2 - The functions 
g n are continuous and g n (t) converges monotonically to x as n tends to infinity for every 
t G [0,T]. Hence g n converges to the constant function r uniformly with respect to t G [0,T]. 
Therefore, for every 1 < k < x, ip^{t) = Y^ = i( < l ) ii'^t4 > k)4 > i converges to T t (pk uniformly with 
respect to t G [0, T]. In particular, the matrix {{^{t), ^y(O)) -f.-i conver g es to I r uniformly 
with respect to t G [0,T]. 

For every n large enough, the vectors ip™(t),. . . ,ip™(t) are linearly independent and can 
be completed to a basis B n {t) = (ip™(t), . . . , ^(t), <p™ +1 (t), . . . , <^"(t)) of span{0i, . . . , n } de- 
pending continuously on t G [0,T]. Let M n (t) be the n x n matrix of the components of 
B n (t) with respect to the basis (0i, . . . , n ). Denote by F n (t) the matrix whose columns are 
the Gram-Schmidt transform of the columns of M n (t). Then, for every n large enough, F n 
satisfies the statement of the proposition. □ 
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We are now ready to prove Proposition 15.11 



Proof of Proposition [3T71 Thanks to Proposition 14. 2 \ it is sufficient to prove that, for every r 
in N, e > 0, and every continuous curve T : [0, T] — > U(H), there exist T u > 0, a continuous 
increasing bijection s : [0, T] — > [0,T U ], and a piecewise constant control function u : [0,T„] — > 
[1/5, +oo) such that the propagator T" of system (14. ip satisfies 

1 1 {^f^}!- 1(^x^0) 1 1 < £ , 

for every (p £ span{0 1; . . . , <p z } with ||0|| = 1 and every t £ [0, T], j £ N. 

By Proposition 15.71 there exist n > r and a continuous curve F n : [0, T] — >■ SU(n) such 
that |(0 j7 f t fc ) - (0j n) , F n (t)0j; n) )| < e/3 for every t in [0,T], 1 < k < r, and j £ N. 

By Proposition 15.4} there exists an admissible trajectory x : [0,Xy — > SU(n) of system 
(E n ) with initial condition I n satisfying 

\\x(s(t)) - F n (t)\\ < \ for every t £ [0,T]. 

Finally, by Proposition 15. 6[ there exists a piecewise constant function u : [0,Te] — > 
[1/5, +oo) such that the propagator T" of (14. ip satisfies 

||(4 n) ,^)n n 0)|-|(^,Tr(0)>||<| 

for every <fi £ span{0 1; . . . , <p z } with ||0|| = 1 and every t £ [0, Tg], j £ N. □ 
5.6 Estimates of the L 1 norm of the control 

We derive now estimates of the minimal L 1 norm of the control u whose existence is asserted 
in Proposition 15.11 We focus here on the physically relevant transitions inducing permutations 
between eigenvectors of A. 

The strategy to get L 1 estimates is the following. Recall that, instead of considering the 
control system x = (A + uB)x driven by a piecewise continuous function u : [0,T U ] — > [0,5], 
we have defined the function V(u) : [0, ||it||x,i] — > [1/5, oo) and considered the control system 
x = (V(u)A + B)x. By Propositions 14.11 and 14. 2\ in order to estimate the L 1 norm of u, it 
is enough to estimate the time needed to transfer the system x = (uA + B)x from a given 
source to an e- neighbor hood of a given target. We observe that the time needed to transfer 
x = (uA + B)x from one state to an e-neighborhood of another is smaller than or equal to 
the time needed to transfer system (E n ) between the n-Galerkin approximations of the initial 
and the final condition for n large enough. 

We proceed to the proofs of Proposition 12.81 and Theorem 12.131 

Proof of Proposition \2.8\ . Let S be a non-resonant chain of connectedness for (A, B, [0, 8], <&). 
Choose (j, k) in S and let m in N be such that m > j, k. The solution x : R — > SU{m) of the 
Cauchy problem 

x = u\b jk \ (e^ - e£j°J x, x(0) = I m 

is a trajectory of (S m ) and has the form x(t) = exp {^f^ — e i7''))' matrix 

M := x C&fa) satisfies = ^ if / ^ j, k, M<\>f ] = -<j)^\ M<J)[ m) = <^ m) . In other 
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words, the control system (E m ) can exchange (up to a phase factor) the eigenstates j and k 
of A^ m \ leaving all the others eigenstates invariant, in time „ ,T , . □ 

Proof of Theorem \2.13l Using Proposition 14.51 one may assume that S is a m-connectedness 
chain for every m G N. Let us prove by induction that every permutation of {1, ... , m} is a 
product of at most 2 m ~ 1 — 1 transpositions of the form (j k) with (j, k) in S D {1, . . . ,m} 2 . 
Let fr(n) be the minimal integer such that every permutation of {1, . . . ,n} is the product of 
at most h(n) transpositions of the form (j k) with (j, k) in S m . For every permutation a of 
{1, . . . , n + 1}, either a(n + 1) = n + 1, and a is generated by at most h{n) transpositions, or 
a(n+l) < n+1. In this case, there exists 1 < k < n such that (k, n+1) G S. Since the product 
(k n + l)(/c cr(n + l))cr leaves n + 1 invariant, it is a product of at most h(n) permutations. 
As a conclusion, h(n + 1) < 2h(n) + 1 and since h(2) = 1, we find h(m) < 2 m_1 — 1. 

The time needed for each of the transpositions (j k) with (j, k) in S has been computed in 
Proposition 12.81 The conclusion follows from the estimate v > 2/5 proved in Remark 14.41 □ 

Remark 5.8. The bound given in Theorem 12. 131 does not depend on e. However, it is possible 
that the time T u needed to achieve the transfer of system (12. ip grows to infinity as e tends to 
zero. 



Remark 5.9. Theorem 12.131 could be stated in a more general way. Indeed the result [26| 
Theorem 6.2] gives the existence of a uniform bound on the time needed to steer system (E n ) 
from any linear combination of the first n eigenstates to any other. This fact guarantees the 
existence of a uniform bound for the L 1 -norm of a control steering system ( 12. ip from any linear 
combination of the first n eigenstates to any neighborhood of any other unitarily equivalent 
linear combination of the first n eigenstates. Such time estimates have been given explicitly 
in the case S = N 2 in [TJ Section 5]. This result could be generalized to the case under 
consideration. It is, however, rather technical and involves advanced notions of Lie group 
theory. 

Following the method of [TT| Section 4.5], one can also give a lower bound for the L 1 norm 
of the control. 

Proposition 5.10. Let (A, B,U,&) satisfy (21). For every T in XJ(H), m inN, e > 0, 
8x, . . . , 6 n G R, and every piecewise constant function u : [0, T u ] — > U such that the propagator 
T u of ([21]) satisfies \\e i6k t (0 fc ) - T^(0 fc )|| <e forl<k<m, one has 

„ „ . (<l>j,<t>k) ~ \{(f>j,^(f>k}\ \ -e 
IMIl 1 > sup SU P 
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Notice that while some strong assumptions about the existence of connectedness chains 
are needed in Theorem I2.13[ Proposition 15.101 is valid even if (12. ip is not approximately 
controllable. 



6 Phase tuning 

Based on Proposition 15.11 we shall now complete the proof of Theorem 12.111 proving approxi- 
mate simultaneous controllability (see Proposition 16.11 below) . 
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In order to outline the mechanism of the proof, we treat in a first time the case of a 
single wave function (proving directly the first part of Corollary I2.12p and we then turn, in 
Section \6.2\ to the general case. 

6.1 Phase tuning for the control of a single wave function 

Simultaneous controllability is obtained from Proposition 15. II applied both to (12.1 J) and to its 
time- reversed version. If (A, B, [0, 5], $) satisfies (21) and admits a non-resonant connectedness 
chain, then the same is true for (—A,—B, [0,5],$). Notice, moreover, that, by unitarity of 
the evolution of the Schrodinger equation, if u : [0,T] — > [0,5] steers ipo to a e-neighborhood 
of ipi for the time-reversed control system 

d ±(t) = -{A + u{t)B)iit), u(t)e[0,S], (6.1) 

then u(T — ■) : [0, T] — > [0,5] steers ipi e-close to for the original system (12. ip . 

Take any eigenvector (p k such that X k ^ (its existence clearly follows from the existence 
of a non- resonant connectedness chain) and consider the control u : [0,T] — > [0,5] steering 
ipo to a e-neighborhood of e td cj) k for some 9 G [0,27r). The existence of such a u follows from 
Proposition 15.11 with Y any continuous curve in U("H) from the identity to a unitary operator 
sending ipo into (f) k and t sufficiently large. 

Similarly, there exist u : [0, T] — > [0, S] and 9 G [0, 2n) such that u steers ip\ e-close to e 
for (16. ip . Let r > be such that 

e-V0s) = e*U~ k - 

Hence, the concatenation of u, of the control constantly equal to zero for a time r, and of 
u(T — •), steers ip 2e-close to ip\. 

6.2 Phase tuning for simultaneous control 

Let r be the number of equations that we would like to control simultaneously, as in Defini- 
tion 12.101 The scheme of the argument is similar to the one above. The pivotal role of the 
orbit of {e tA (p] t \ t} is now played by a torus of dimension r. 

The crucial point is to ensure that an orbit of A "fills" the torus densely enough. This is 
formally stated in the proposition below. Recall that a subset Q\ is e-dense in a metric space 
Q 2 if any point of Q 2 is at distance smaller than e from every point of f^. For every m G N 
and kx, . . . , k m G N, define 

T(h, . . . , k m ) = {e e ' A <p kl + ■■■ + e 9mA <p km | 9 U . . . , 9 m G R} 
C(h, ...,k m ) = {e tA <f) kl + ■■■ + e tA <)) km | t G R}. 

Notice that, if ki, . . . ,k m are distinct, T{k\ . . . , k m ) is diffeomorphic to the torus T m . 

Proposition 6.1. Let (A, B, [0, S], $) satisfy (21) and admit a non-resonant connectedness 
chain. Assume that, for every rj > and r in N, there exist r pairwise distinct positive integers 
ki, . . . ,k r such that C(k~i, . . . , k r ) is rj-dense in Tiki, . . . , k r ). Then (12. ip is simultaneously 
approximately controllable. 
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Proof. Take r orthonormal initial conditions ipQ, . . . ,ipQ and r orthonormal final conditions 
. . . Fix a tolerance rj > 0. Take fcj, . . . , k r as in the statement of the proposition. 
According to Proposition 15.11 (with r sufficiently large and T a continuous curve in 111(7/) 
from the identity to a unitary operator sending ipQ to fa for each j = 1, . . . , r), there exists 
a control u(-) steering simultaneously each ip 3 , for j = 1, . . . ,r, //-close to e ej fa. for some 
0i,..., 9 r G R. Similarly, applying Proposition 15.11 to the triple (—A,—B,[0,5],<&), there 
exists a control u : [0,T] — > [0,5] steering simultaneously ip[ ?7-close to e * fa for ( 16. ip . for 
j = 1, . . . ,r and for some 0\ , . . . , r G R. 

Since the positive orbit of A passing through 5^/=i e° iA fa. is 77-close to X^=i e6iA ^kj ^ 
T(ki, . . . , k r ), then the concatenation of «(•), a control constantly equal to zero on a time 
interval of suitable length, and u(T — ■) steers each ip J 3?7-close to ip[ for j = 1, . . . ,r. □ 

Remark 6.2. As it follows from the proof above, in Proposition 16.11 the hypothesis of exis- 
tence of a non-resonant chain of connectedness can be replaced by the weaker hypothesis that 
both (12. ip and (16. ip are simultaneously controllable up to phases in the sense of Proposi- 
tion ECU 

We are left to prove the following. 

Lemma 6.3. // A has infinitely many distinct eigenvalues then for every 77 > and r G 
N, there exist r distinct positive integers ki,...,k r such that C(ki, . . . , k r ) is rj-dense in 
T{k\, ■ ■ ■ , kr) . 

We split the proof of Lemma 16.31 in four cases. 

6.2.1 dimQ (span Q (A fc ) fceN ) = 00 

This case is trivial: it is enough to take fci, . . . , k r in such a way that A^, . . . , \j. r are Q-linearly 
independent. Then C(ki, . . . , k r ) is dense (hence, 7/-dense for every 77 > 0) in T{k\, ■ ■ ■ , k r ). 

6.2.2 The spectrum of A is unbounded 

The most physically relevant case is the one in which the sequence (\k)k&~N is unbounded. 

Fix k± such that A^ 7^ 0. Take then ^2 such that |A^ 2 | S> |A^J in such a way that the 
orbit C(k\, fc 2 ) is //-dense in T{k\, k-z). By recurrence, taking |A^ m | 3> |Afc m J for m = 2, . . . ,r 
we have that C{k\, . . . , k m ) is 7/-dense in T{ki, . . . , k m ). Indeed, by the recurrence hypothesis, 
C(ki, . . . , k m _i) + T(k m ) is 77-dense in T(k~i, ■ ■ ■ , k m ) = T{k\, . . . , k m ~i) + T(k m ) and the choice 
of Ag m is such that C(k\, . . . , k m ) is 77'-dense in C(ki, . . . , fc m -i) + T(k m ) for rj' arbitrarily small. 

6.2.3 The spectrum of A is bounded and dimQ (spanq(Afc)fcgN) = 1 

The existence of a connectedness chain implies that there exist infinitely many pairwise distinct 
gaps between eigenvalues of A. Hence, the set of eigenvalues of A has infinite cardinality. 

Since all the eigenvalues of A are Q-linearly dependent, we may assume, without loss of 
generality, that Xj is rational for every j G N. Up to removing the eigenvalues equal to zero, 
if they exist, and changing the sign of some eigenvalues we can also assume that 

Aj = ~r~i gcd(a i ,6 j ) = 1, aj,bj > 
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for every j G N. The boundedness and the infinite cardinality of the spectrum of A imply 
that the sequence (bj)j G j$ is unbounded. 

In order to prove Lemma 16. 3[ let us make some preliminary considerations. For every 
m G N and fci, . . . , k m G N, denote by r(ki, . . . , k m ) the minimum of all t > such that t\ k 
belongs to N for every j = 1, . . . , m. Equivalently said, 27rr(fci, • • • , k m ) is the period of the 
curve s i— > e sA (j) kl + ■ ■ ■ + e sA( Pk m - Notice that, if A^, . . . , X km are integers, then r(fci, . . . , k m ) = 
l/gcd(a fcl , . . . , a km ). In general, 

T{k 1 ,...,k m ) = —- f — J^ m r . (6.2) 

gcd {a kl U j=1 b k .Uf =l+1 b kj ) 1 ^ m 

The following lemma guarantees that, for any choice of k±, . . . , fc m _i, we can select k m in 
such a way that r(ki, . . . , k m ) 3> r(fci, . . . , fc m _i). 

Lemma 6.4. Let m G N. For every fci, . . . , k m G N, there exists c = c(k\, . . . , k m ) > such 
that, for every k m+ i G N ? 



r(fci,...,fc m+ i) ^ b 
t(/ci, . . . , /c m ) 



^ 1 



Proof. From equation taum, 



r(fc 1 , . . . , fc m+1 ) , gcd (a fe! n;. = 1 1 6 fe .n^ +1 6 fc .) 1 



<Km 



r(/c 1; ...,k m ) Km+1 gcd (a^ny^ng^^) 



KKm+l 



The proof consists, then, in showing that gcd { a h^ l jJibk j ^ 7 j l =i+^ :>k j)i<i<m+i * s bounded 
from above by a constant independent of k m +i- First, notice that 

gcd Knfc^^n^iA,)^^ < gcd (d^,^^) =: r. 

Set Ci = a k JY^ =2 b kj and C2 = n^o^. and notice that they do not depend on k m+ i. 
Write T as 

7i/3 m+ i = T = 7 2 a m+ i 

where 71 and 72 divide Ci and C2, respectively, while a m +i and /3 m+ i divide a km+1 and b km+1 , 
respectively. Since a km+1 and Ofc m+i are relatively prime, then the same is true for a m+ \ and 
f3 m+ i. Therefore, a m+ i divides 71. Hence, T = a m +i72 < 7172 < □ 

We show now how to choose kx, . . . ,k r as in the statement of Lemma 16.31 We proceed 
by induction on r. The case r = 1 has already been treated in Section 16.11 Assume that 
k±, . . . , fc r _i are such that C(k\, . . . , k r -i) is ?7-dense in T(ki, ■ ■ ■ , k r -i)- Hence, for every choice 
of k r , the set C(ki, . . . , fc r _i) + T(k r ) is ?7-dense in T(ki, . . . , £v). 

We are left to show that, for a suitable choice of k r G N\{fci, . . . , fc r _i}, the set C(ki, . . . , k r ) 
is r/'-dense in C(ki, . . . , fc r _i) + T(k r ) for 77' arbitrarily small. Recall that C(k\, . . . , fc m ) is 
the support of the curve s 1— > e sA (p kl + • • • + e sA (f> krn , whose period equals 27rr(fci, . . . , k m ). 
Therefore, C(k\, . . . , k r -i) is the projection of C(ki, . . . , k r ) along <p kr on 0^. Hence, for 
every ip G C(ki, . . . , k r -i) the cardinality of the set (ip + T(k r )) H C(&i, . . . , fc r ) is equal to 
r(ki, . . . , k r )/r(ki, . . . , fc r _i). In particular, + T(k r )) PI C(fci, . . . , fc r ), which is regularly 
distributed, is 27ir(ki, . . . , fc r _i)/r(fci, . . . , fc r )-dense in ^ + T(k r ). 

Lemma [6.41 and the unboundedness of the sequence (o^OjeN allow to conclude. 
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6.2.4 The spectrum of A is bounded and 1 < dimg (spang (Afc)fc e N) < oo 

Let m = diniq (spanq(Afc)fc G N) and fix a Q-basis . . . ,/i m of spanq(Afc)fe e N- 
Let 

m 

Xj = ^2 a i^i a ], . . . ,af e Q. 

i=i 

There exists I G {1, . . . , m} such that the cardinality of {a*- | j G N} is infinite. (Otherwise, 
the set of eigenvalues of A would be finite.) Without loss of generality, 1 = 1. 

The results of Sections l6.2.2l and l6.2.3l implv that, for every rj > 0, there exist k±, . . . , k r G N 
such that {(e** Mia, =i , . . . , e itMia *' ) | t G R} is //-dense in T r . We are going to show that 
C(ki, . . . , k r ) is rmry-dense in Tiki, . . . , k r ) or, equivalently, that {(e ltXk ^, . . . , e ttx ~ k r) \ t G R} 
is rm?7-dense in T r . 

Up to a reparametrization, we can assume that al G Z for every j = 1, . . . ,r and / = 
l,...,m. 

Fix (e %ei , . . . , e j6,r ) in T r . The choice of ki, . . . , k r G N guarantees the existence of t G R 

such that ||(e* <Ml0fc i, . . . , e 1 ^ 1 "^) — (e 101 , . . . , e j6lr )|| < 77. Because of the Q-linear independence 
of Hi, . . . , /i m , there exists < G R such that 

IKe^.e*** • • • , e^») - (e*" 1 , 1, . . . , 1)|| < V — -. 

max{|ay | j = 1, . . . ,r, / = 1, . . . ,m| 

it(j,ia l - 

In particular |e ^ — 1| < 77 tor every I = 2, . . . , m and every j — 1, . . . , r. 
Hence 

|| (e^i , . . . , ) - (e iei , . . . , ) || < || (e ltx h ) - (^ ia h e if " ia i ) || 



ie €/xia i)-(e^,...,e ie ^ 

it/iiaj_ iifimCtji' itfiictj; 



< > I e i ■ ■ ■ e i — e j \ + rj 

i=i 

< (r(m — 1) + 1)77. 

This concludes the proof of Lemma 16.31 and of Theorem 12.111 □ 



7 Example: Infinite potential well 

We consider now the case of a particle confined in (—1/2, 1/2). This model has been extensively 
studied by several authors in the last few years and was the first quantum system for which 
a positive controllability result has been obtained. Beauchard proved exact controllability in 
some dense subsets of L 2 using Coron's return method (see [51 E] for a precise statement). 
Nersesyan obtained approximate controllability results using Lyapunov techniques. In the 
following, we extend these controllability results to simultaneous controllability and provide 
some estimates of the L 1 norm of controls achieving the transfer between two density matrices. 
The Schrodinger equation writes 
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with the boundary conditions 1/2, t) = ip(l/2,t) = for every t G R. 

In this case "H = L 2 ((—1/2, 1/2), C) endowed with the Hermitian product (^1,^2) — 
I-1/2 1 Pi{, x )' l l J 2{.x)dx. The operators A and I? are defined by Aip = for every ip in 

D(A) = (H 2 n tfo 1 ) ((-1/2, 1/2), C), and 5^ = zx<0- 

Unfortunately, due to the numerous resonances, we are not able to apply directly our 
results to system (17.11) . A classical approach is in this case to use perturbation theory. Indeed, 
consider for every 77 in [0, 5], the operator A v = A + t]B. The controllability of system (17.11) 
with control in [0, 5} is equivalent to the controllability of 

dib 

-j- = A^ + vBijj 
at 

with controls v taking values in [— rj, 5 — 77]. Since the perturbation 77 i-> A rj is analytic, the 
self-adjoint operator A v admits a complete set of eigenvectors ((f>k{'r]))k£'N associated with the 
eigenvalues (iXk{v))k£N, with 0*. and analytic (see |T7]). For 77 = 0, 



= (0) 



x 1 — y \[2 cos(kirx) when k is odd 
x H- \/2 sm(k-Kx) when k is even 



is a complete set of eigenvectors of A associated with the eigenvalues i\k{0) = —i^-. Fol- 
lowing [2 Proposition 2.3], one can compute the 2-jet at zero of the analytic functions 

as 77 goes to zero. For every k\, ^2,^1,^2 in N, since \ is transcendental on Q, 



^2 Pi P2 



1 1 . . 1 _ 1 

P'2 



73 — tk — -z — -r 

K l rc 2 




1 1 w 1 

^ — Z2" I 3 T 3 



fcl =pi 

k 2 = P2- 



k Pl P2/ \PT P2, 

Hence, the 2-jets of the gaps between two eigenvalues are all different at zero. Recall that 
B(f>k +1 ) ^ 0. An argument similar to the one in [HI Proposition 6.2] ensures that for 
every e > there exists < 77 < e such that {(j, j + 1), (j + l,j) | j G N} is a non-resonant 
connectedness chain for (A v , B, [—77, <5 — 77], $). 

We can now apply Theorem 12.111 to obtain simultaneous approximate controllability for 
(17. ip and Theorem 12. 131 to get L 1 estimates. For instance, for every s > there exist T > 
and a piecewise constant function u : [0,T] — > [0,5] such that the propagator at time T of 
(17. ip exchanges (up to a correction of size e) the density matrices 

^>i(0)<M0)* + ^ 2 (O)0 2 (O)* and ^ 2 (O)0 2 (O)* + ^i(O)0i(O)* 

and 

7T 9-7T 3 

w L i < — n — r~, rrrr = ~ 20.2656 . 

11 l|L " 2^1(^(0), B<j> 2 (0))\ 32v 
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The control time T satisifies T > | 
satisfy 



u\\li and, from Proposition 15.101 the control u has to 



^Hl 1 > (1 — ^) max 



2\/37r 



(l-e)5.5323. 



V2tt 2 -3 



8 Orientation of a bipolar molecule in the plane by means 
of two external fields 

We present here an example of approximately controllable quantum system. (See [2"Tl ESI 
and references therein.) It provides a simple model for the control by two electric fields of the 
rotation of a bipolar rigid molecule confined to a plane (see Figure [2]). Molecular orientation 
and alignment are well-established topics in the quantum control of molecular dynamics both 
from the experimental and theoretical point of view. 

The model we aim to consider can be represented by a Schrodinger equation on the circle 
S 1 = R/27rZ, so that "H = -^ 2 (S 1 , C). In this case A = i d 2 /d8 2 has discrete spectrum and its 
eigenvectors are trigonometric functions. The controlled Schrodinger equation is 



We assume that both u\ and U2 are piecewise constant and take values in [0, 5], 5 > 0. 

Notice that system (18. ip is not controllable if we fix one control to zero. Indeed, by parity 
reasons, the potential cos(#) does not couple an odd wave function with an even one and the 
potential sin(6') does not couple wave functions with the same parity. 

For every a 6 S 1 , let us split H as © H", where H° (respectively, is the closed 
subspace of H of even (respectively, odd) functions with respect to a. Notice that and 

are Hilbert spaces. A complete ortho normal system for Ji® (respectively, H") is given 
by {cos(fc(- — tt))/v / ^'}felo (respectively, {sin(fc(- — &)) / y/n}kLi) • Let us write ip = ip" + ip°, 
where G U a e and ip™ G H". 

Our first result states that, once a is fixed, the even or the odd part of a wave function ijj 
can be approximately controlled, under the constraint that their L 2 -norms are preserved. 

Lemma 8.1. Let ip G "H and a G [0, 7r/2]. Then for every e > and every ip G "H" such that 
W^elW — W^Wh, there exists a piecewise constant control u = (wi, 1 ^) steering xjj to a wave 
function whose even part with respect to a lies in an e -neighborhood ofip. Similarly, for every 
e > and every ip G "H" such that ||% = H'i/'llw) there exists a piecewise constant control u 
steering %p to a wave function whose odd part with respect to a lies in an e -neighborhood ofip. 

Proof. The idea is to apply Theorem 12.61 to a subsystem of ( 18. ip corresponding to the choice 
of a subclass of wave functions and a subclass of admissible controls. 
More precisely, let 

W a = {u : R — > [0, 5} 2 piecewise constant | u(t) is proportional to (cos a, sin a) for all t G R}. 




(8.1) 



26 



The control system whose dynamics are described by (18.11) with admissible control functions 
restricted to U n can be rewritten as 



' dt 



d 2 



do* 



+ v(t) cos(9 — a 



M ip(9,t), v G (o,5y/l + minjtan a, cotan a}' 



Notice that the spaces "H" and are invariant for the evolution of ( 18. 2p 

d 2 



(8.2) 

whatever the 



choice of v — v(-), since they are invariant both for A = i-^p and for the multiplicative 
operator B a :U^U defined by (B a (p)(8) = -icos(6 - a)(p{9). 

We shall consider ( 18. 2 p as a control system defined on (the second part of the statement 
of the lemma can be proved similarly by considering its dynamics restricted to H"). 

Denote by A" and B° the restrictions of A and B a to H° . Choose as orthonormal basis 
of eigenfunctions for A" the sequence defined by (f)k(0) = cos(k(6 — a))/-^ for k G N. The 
eigenvalue of A associated with (p^ is iX^ = —ik 2 . 

Then (<pj, B"(pk) ^ if and only if \k — j\ — 1. We take as connectedness chain the set 
{(k,j) G N 2 | \k — j\ = 1}. Since X^+i — = —2k — 1, then the connectedness chain is 
non-resonant. Theorem 12.61 implies that ( 18. 2 p can be steered from ip® to an e-neighborhood 
of any ip G such that HV^II = H^H by an admissible control v(-). The conclusion follows 
by applying u(-) = (v (•) cos «,?;(•) sin(a)) to ( 18. ip (since and are invariant by the flow 
generated by «(•))• □ 



The main result of this section states that (18. ip is approximately controllable. 
Proposition 8.2. System ( 18. ip is approximately controllable. 

Proof. It is enough to prove that every wave function of norm one can be steered arbitrarily 
close to the constant \ j\phx. Indeed, if ip and ip have norm one, if the control u(-) steers the 
initial condition ip e-close to the constant 1/a/27t, and if «(■) steers the conjugate of ip e-close 
to the same constant 1 / V2%, then the concatenation of u and of the time reversed of u steers 
ip 2e-close to ip. 

Fix ip G "H of norm one, a tolerance e > 0, and choose a G (0, 7r/2). Fix s = e(l — l/y/2). 
Then, according to Lemma IHTTj ip can be steered to a wave function ip such that \\ip — ipi\\ < e, 
where ip\ is of the form 

m 



2tt 



If ||0i|| is smaller than e/2 then we are done. Indeed 
which implies that 1 — \\ip" \\ < e 2 /4. Then 



+ <p u with fa G K- 

i has L 2 -norm equal to 



Set 112 



iP- 



'2tt 



< 





1 


Tpl- 




V2^ 



2^/l-\\iP«\\<e+-j=e = e. 



<e+y/(\m\-l) 2 + 

1 



Assume then that > e/2 and consider, for every /3 G S 1 , 

characterize Tg in terms of the coefficients of the representation 



i)f || 2 . We can 



k=l 



sm(k( ■ — a)) 
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Indeed, 



4 = E I ^i('), c os(fc(- - /3))/V5F) | 2 = E IK sin(fc( • - a))/V5F, cos(fc(. - /3))/v^) | 

k=l k=l 

oo 

|a fc sin((/3 - a)fc)|' 



k=i 

There exists c > independent of k and a such that 

sin 2 ((/3 - a)jfe)d/3 > c, 



o 



hence, 

oo 



/ 2 ^|a fc | 2 sin 2 ((/3-a)£;)d/3>cJ] 

"0 t— i i— 1 



a k \ 2 = c||0i|| 2 , 



k=l k=l 

from which we conclude that there exists (3 G (0, 7r/2) such that 

^>-||0iH 2 >f^ 2 . 
Notice now that the even part of ^ with respect to /3 has norm 



2 



i)i\\ = y/\m\ 2 +r P >^\m\>+^. 

Repeating the same argument as above replacing ip by ipi we conclude that it is then possible 
to steer ip at a distance smaller than e from the sum ip2 of a positive constant function of norm 



larger than y HV'e || 2 + and a function 02 6 %f of norm ||02|| < ||0i||- If ||02|| < s/2 then 
we are done. Otherwise, since the improvement in the size of the constant is bounded from 
below by a quantity that does not depend on <pi, we can iterate the procedure finitely many 
times up to guaranteeing that the final wave function is £-close to the constant 1 / y/2n. □ 

A Appendix: relations between controllability notions 

Proposition A.l. Let (A,B,U,$) satisfy (21). Then 

(i) for every r in N, e > 0, and T in U(%) there exists a piecewise constant control 
u : [0,T] — y U such that such that ||T0fc — ^r^kW < £ for every 1 < k < r, 

implies 

(ii) (12.1 p is controllable in the sense of densities matrices 
which implies 

(Hi) for every r in N, £ > 0, and T in U("H) there exist 9i,...,9 r G R and a piecewise 
constant control u : [0, T\ — > U such that \\e l6k T((pk) — Y^(0fc)|| < £ f or every 1 < k < r. 
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Proof. Let us prove that (i) implies (ii). Fix two unitarily equivalent density matrices p 
and p\. Write p = YlkeN Pk v k v k with (^fc)fceN an orthonormal sequence in H. and (Pk)keN a 
sequence in ^ 1 ([0, 1]) such that Pj. = 1. By assumption, there exists T in U("H) such that 

Pl = fpot* = E fceN ^tK)tK)*. 

Let m in N be such that ^2 k>rn Pt < £ and r in N be such that \\vj — Ylk=i v o) < £ 
for every j — 1, . . . ,m. By hypothesis, there exists a piecewise constant function u : [0, T] — > U 
such that ||T(0fc) — Ty(0k)|| < e/r for every 1 < k < r. Hence, for j — 1, . . . ,m, 



\T(v 3 )-T u T (v 3 ] 



£<^,«i>(f(^)-T^(^)) 



fc=i 



< 



£|| (0 fc ,^)(f(0 fc )-T^(0 fc ) 



fc=i 



+ 



K k=r+l 



<e + 2 



k=l 



< 3e. 



Then, recalling that for every 9 in R, a, b in "H, \\aa* — bb*\\ < \\aa* — e t6 ba* + e ld ba* — bb*\\ < 
||a||||a-e ie 6|| + \\b\\\\e w a* - b* \\ < (\\a\\ + \\b\\)\\a - e te b\\, we get 

m 

\mpoTf - pi\\ < ^m^ j )r^(v j y-f(v j )f( Vj y\\ + 2e 



3=1 
m 



< E^(H t t(^)II + 11^)11)11^-) - t( Vj ] 

3=1 

< 6e + 2e = 8e, 



2e 



which concludes the first part of the proof. 

Assume now that (ii) holds true. Fix e, r G N, and T as in the hypotheses. Choose 
eti, . . . , a r G R such that < a\ < a 2 < ■ ■ ■ < a r and Y7k=i a k = !• Define the two unitarily 
equivalent density matrices 



k=l 



ak<Pk<Pk 



and Pi = y^Qfcf (0 fc )f ( 



k=l 



By assumption, there exists a piecewise constant it : [0,T] — )■ [/ such that 

||pi-T£p T£l <Ce, 



(A.1) 



where C = min{aj, |afe — aj| | I < j,k,l < r,k ^ 0/2. Choose 1 < A; < r and test flA.ll) on 



aJ^ -^a fc (TW,T^ )t( 



fc=i 



< Ce. 
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Since (T((j) k )) k is an Hilbert basis of H, then 



r u T( p ko = J2^(^),r u T 4> k0 )f( 



k=l 



and hence 



J>*o - a k )(f((f> k ), T^ k0 )f (4> k ) + a k0 J](f (cf> k ), T^ feo )T (<f> k ) 



k=l 

In particular we have 



k>7 



< Ce. 



^<f(0 fc ),T^ feo }T(0* 

fc^fco 



£ 

< 2" 



For £ small enough, this leads to |(T(0fc o ), Y^0fc o )| > \J\ — e 2 /A > 1 —e/2. Hence, there exists 
6 ko in R such that ||e**of (0 fc() ) - T£(0 fco )ll < e. □ 



Proposition A. 2. Under the hypotheses of Proposition I7Q1 and if, moreover, (<frj,B<j) k ) is 
purely imaginary for every j,k £ N, i/ien (iii) implies (ii). 

Proof. First of all, because of the hypotheses on B, system (12.11) satisfies the following time- 
reversibility property: If ip : [0, T] — > 7-L is a solution of (12. ip with control u(-), then <p(t) = 
SfcLi ^{T — t))(f) k is a solution of ( 12. ip with control u(T — •). 

Fix two unitarily equivalent density matrices po an d pi- Write po = SfceN PkVkV^ with 
(ffc)fceN & n orthonormal sequence in "H and (Pfc)fcgN a sequence in £ l ([0, 1]) such that J^ fcgN P k 
I. Let f in U(ft) be such that Pl = Tp T* = £ fceN P fc f (u fc )f 

Let r £ N be such that Ylk>r < £ - By hypothesis there exist a control u : [0, T u ] — > [0, 5] 
and 9i,...,9 r £ R such that ||e i6,fc X^li ( < A?>' y fc) ( A? ~ ^T„(^fc)ll < e f° r ever Y 1 < & < r. By 
time reversibility the control u(T u — •) steers e~ l9k v k e-close to cf) k . Now, there exist a control 
w : [0, T w ] y [0, 5] and 9i, . . . ,9 r such that 

||e^f(^)-T^(0 fc )||. 

Let £t be the concatenation of the controls u(T u — •) and w.Then 

\\e^?(v k )-r%(e-^v k )\\<2E. 

Finally, 



\\n u po?% - Pill < J2 p ^ T U^ i9kv k) T U e ' i6kv kT - f(v k )f(v k y\\+ 

k=l 

oo 

+ E Pk\\rl(v k )r^(v k y -t(v k )t(v k y\\ 

k=r+l 
r 

< 4e Pk + 2e < 6e. 

This concludes the proof. 



k=l 



□ 
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